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The occurrence of a bubble, due to an inversion of Si/2 state with the state usually located above, 
is investigated. Proton bubbles in neutron-rich Argon isotopes are optimal candidates. Pairing 
effects which can play against the bubble formation are evaluated. They cannot prevent bubble 
formation in very neutron-rich argon isotopes such as ®®Ar. This pleads for a measurement of the 
charge density of neutron-rich argon isotopes in the forthcoming years, with the advent of electron 
scattering experiments in next generation exotic beam facilities such as FAIR or RIBF. 
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I. INTRODUCTION 

The occurrence of bubble nuclei has been studied by 
Wilson in 1946 1], in order to describe low- lying excita- 
tions as oscillations of nuclear bubbles. The invoked rea- 
son for the bubble existence was an overestimated "sat- 
uration of the nuclear forces" . Bethe and Siemens have 
also investigated this possibility in the sixties j2|- The 
stability of the bubble was studied using the liquid drop 
model and this analysis showed that bubbles were not 
leading to the minimum energy configuration, compared 
to prolate shape with uniform density. During the sev- 
enties and eighties, bubble nuclei were also investigated, 
still with liquid drop based models 0] or Thomas-Fermi 
approaches More recently bubbles were predicted in 
hyperheavy nuclei Q with Z > 120. In this case, self- 
consistent Hartree-Fock-Bogoliubov (HFB) calculations 
have shown that the ground state corresponds to a bub- 
ble configuration. This effect is due to an interplay be- 
tween the Coulomb interaction and the nucleon-nucleon 
interaction in those very heavy systems. The optimum 
spatial position of the bubble has also been investigated 
0. However the existence of such heavy nuclei is specu- 
lative, and bubbles have not been detected yet. It should 
be noted that bubble effects may be a general feature 
in nuclei: they have also been predicted in superheavy 
nuclei around Z=120 using mean-field calculations. 

A qualitatively different interpretation for the bubble 
existence is invoking a specific microscopic quantum me- 
chanical effect: the s-states being the only wave functions 
with non zero value at r=0, the depopulation of this level 
leads to a depletion of the density in the center of the 
nucleus. Such an interior-peaked shape of the 3s wave 
function has been measured using electron scattering on 
2"6pb and 205T1 |8|]. However, the depletion in the inte- 
rior of these nuclei remains small since in heavy nuclei the 
2si/2 state is already filled and contributes to the den- 
sity at r=0. In an early Hartree-Fock-|-BCS calculation a 
strong depletion in the interior of "^^Ar and ^"''Hg nuclei 
was predicted @, which the authors defined as a bubble 
occurrence. However, their effective nucleon-nucleon in- 
teraction was not considered as sound, and it was shown 
that using more elaborated interactions prevent bubble 
formation in '^^Ar fTol|. It should be noted that the key 



point for the bubble formation in the s-d region is the in- 
version of the proton 2si/2 and ld3/2 states with respect 
to the standard shell model. 

The emergence of experimental studies on exotic nuclei 
gathered numerous signals of shell structure modification 
these last decades (see e.g. It is therefore le- 

gitimate to suspect that in more neutron-rich nuclei, the 
s and d state crossing may occur. Such an inversion was 
recently predicted in the ^^Ar nucleus using relativistic 
mean field (RMF) calculations. A strong depletion in the 
center of the proton density was obtained [ij] . However 
pairing effects could preclude the bubble effect due to the 
occupancy of 2si/2 state. Such a study has been extended 
to N=28 isotones The mechanism of this inversion 
in neutron-rich nuclei has been analyzed in Ref. [l6l| , and 
it is due to the conjunction of several factors: the spin- 
orbit potential modification, the tensor force, as well as 
the lowering of the proton potential and the extension of 
the proton density due to the presence of a neutron skin. 
Predictions for more neutron-rich nuclei show an increase 
of the gap between the inverted s and d states, strongly 
supporting the formation of proton bubbles in ^""^^Ar. 

In this paper we analyze the general features of bubble 
occurrence in exotic nuclei, on an illustrative case. The 
bubble occurrence is a direct consequence of the Si/2 de- 
population, which could be obtained either by lowering 
the proton number, or by a level inversion involving an 
s state. This work is devoted to the latter case. The 
possible candidates are reviewed in Section II, showing 
that proton bubbles in very neutron-rich Argon isotopes 
are the most likely. Several experimental methods to de- 
tect such a bubble are then investigated : charge density 
measurements (Section III) and direct reactions (Section 
IV). 



II. CANDIDATE NUCLEI 

A. Possible candidates 

Bubble candidates are nuclei where the s state is de- 
populated. In some cases this occurs due to an inversion 
between the s state and the one usually located above. 
Therefore, either proton or neutron states inversion be- 
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FIG. 1: Proton densities of *^Ar calculated with the SkI5 
interaction (full lines) and SLy4 [23| (dashed lines) in the HF 
approach. The ^^Ar proton density calculated with the SkI5 
interaction is also shown (dashed-dotted lines) 



tween (2si/2,ld3/2) states or (3si/2,lliii/2) states are pos- 
sible candidates. In the proton bubble case, these inver- 
sions correspond to Ar and Hg isotopes. In the neutron 
bubble case, they correspond to the N=18 and N=80 iso- 
tones. A first experimental signal for an inversion may 
be given by the J'^ of the ground state of the odd nuclei 
having one additional nucleon compared to the bubble 
nuclei (N or Z equals to 19 or 81). A J''=l/2+ together 
with a spectroscopic factor close to one is an indication 
for an Si/2 state located above the usual ld3/2 or lhj^j^/2 
state. In the case of the Potassium isotopes, such val- 
ues have been measured in ''^K [13], showing that the s-d 
states are almost degenerated. An (2si /9,ld3/2) inversion 
is even predicted by RMF calculations [13]. In Ref. p^ . 
the study along the Calcium isotopic chain shows that 
the s-d states inversion is strengthened for very neutron- 
rich nuclei (A>58). Hence the study of very neutron-rich 
Argon isotopes such as ^*Ar using the HF method pro- 
vides an illustrative study of the bubble manifestation. 
Fig. [T] shows the proton densities of ^*Ar, calculated 
with spherical HF model using the SkI5 [isll force. The 
numerical details are the same as in Ref. [la|. A strong 
depletion in obtained in the center of the nucleus, with 
this SkI5 parametrization. It should be noted that the 
large s-d state inversion for very neutron-rich nuclei in 
the ^*Ar region is also obtained with several other inter- 
actions [31 as well as shell-model calculations [Toj . 

In the case of the N=18 isotones, the measured J'^ 
values of the ground states of the N=19 isotones are 
J'^=3/2'^ [2l[, showing no clear signal for an inversion. 



On the contrary all the measured J'^ of the Z=81 iso- 
topes are 1/2+ in the ground state, showing a systematic 
inversion leading to a possible proton bubble in the Hg 
isotopes. However, as stated above, the depletion in the 
center of the proton density of the Hg isotopes is less 
pronounced than for the Argon isotopes, since in heavy 
nuclei the 2si/2 state is already filled and contributes to 
the density at r=0. For the N=81 isotones, J'^ are mea- 
sured for only 3 nuclei and J^=l/2+ is found in ^"'^Gd 



and ^ Dy [2lJ. In this case also, the bubble effect m 
heavy nuclei is expected to be less important than in in- 
termediate mass nuclei. To summarize, the optimal man- 
ifestation of a bubble could be for protons in neutron- rich 
Argon isotopes. We will therefore focus on these nuclei 
in the following. 

It should be noted that deformation can play against 
a bubble occurrence: the degeneracy removal attenuates 
the bubble effect generated by the inoccupancy of the 
2si/2 state. In the case of '^^Ar several N=28 studies used 
deformed mean-field models [12, IH, [l^], and also the 
Generator Coordinate Method (GCM) based on Gogny 
force [1^. In the HFB calculations, ^^Ar is predicted as 
a soft nucleus, either spherical or with a small deforma- 
tion parameter in the ground state. More neutron-rich 
Argon isotopes are also predicted spherical [13]. In the 
GCM approach, the ground state of ''^Ar is predicted as 
a possible shape coexistence state. In the following we 
will consider the spherical case for both "^^Ar and ^^Ar 
nuclei in order to study the bubble hypothesis in this 
framework. 



B. Pairing effects 

Pairing effects may hinder the bubble formation since 
scattered pairs could populate the Si/2 state, decreasing 
the depletion in the center of the nucleus. The occu- 
pancy factor of the s-state due to pairing effects may 
be modified. The study of the proton pairing effect in 
Argon isotopes is a delicate task : it is occurring in a 
two- hole state from a doubly magic nucleus. It is known 
that the HFB approximation is not well designed to such 
a situation, leading to an overestimation of the pairing 
effects. For instance, we have performed HFB calcula- 
tions on 46,68^j. ygjjjg SkI5 interaction and a delta 
density-dependent pairing interaction, of the form : 



Vpair — Vq 



V 



Po 



S (ri - r2) 



(1) 



with ?/=l, a=Q.5 and po=0.16 fm''. The magnitude 
Vo=-330 MeV.fm'^ of the pairing interaction is obtained 
by reproducing the two-proton separation energy of Ar. 
Fig. 2 displays the proton density of ^^Ar, showing a 
reduction of the bubble effect, but still a depletion is 
predicted in the center of the nucleus. 

In *^Ar the pairing interaction plays an even more dra- 
matic role: the s state occupation is predicted to be 54%. 
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FIG. 2; Proton densities of '^''Ar (dashed lines) and ''^Ar 
calculated with the SkI5 interaction in the HFB approach. 



The corresponding proton densities, with and without 
pairing effects, are displayed on Fig. 1 and 2, respec- 
tively. As expected the ld3/2 state is depopulated, due to 
the pairing effect. Shell- model calculations also predict 
around 50% occupation probability for the 2si/2 state 
[T^ . However, it should be noted that the 2si/2 and 
ld3/2 occupation probabilities have not been measured 
yet, leaving a small possibility for a bubble occurrence in 
46 Ar. 

In Ref. the splitting of the neutron p-states in 

''^Ar is weakened by the proton density depletion in the 
center of the nucleus, which modifies the spin-orbit term 
in the interior: the proton depletion is related to the neu- 
tron splitting of the p orbits measured in Ref. [l^, [2^ . 
The weakening of the neutron p states splitting with the 
occupancy of the 2sl/2 orbital has been recently studied 
in [l^ from Si to Ca N=28 isotones. A proton deple- 
tion is predicted for nuclei having a low 2si/2 occupation 
probability. 

The presence of a proton bubble in ''^Ar is subject to 
significant uncertainties, due to the dependence on the 
pairing and the Skyrme interactions. Hence, experimen- 
tal indications of a bubble in this nucleus are of relevant 
importance. The probability of a proton bubble mani- 
festation increases when going from ^^Ar to much more 
neutron-rich nuclei such as 6*^"^* Ar. As mentioned above, 
several interactions used in Ref. predict an inversion 
between the 2si/2 and ld3/2 proton states, with a large 
gap around ''°Ca, showing that the proton bubble pre- 
diction for very neutron-rich Argon isotopes is a general 
feature. 



In the following we choose ^^Ar as an illustrative ex- 
ample. We will therefore consider the bubble hypothesis 
through the HF predictions in ^^Ar using the SkI5 in- 
teraction (Fig. [1]). The no-bubble hypothesis can be 
obtained either by SLy4-HF (Fig. [J) or by SkI5-HFB 
(Fig. [2]) calculations. We have checked that the two 
latter cases give very similar results: in Ref. 16^ HF cal- 
culations using SLy4 do not predict any inversion of the s 
and d states, implying no bubble. Hence we will describe 
the no-bubble hypothesis using the HF calculations with 
the SLy4 interaction. The case of ''^Ar is chosen only for 
illustrative purposes since experimental data are avail- 
able on this nucleus. All conclusions below remain valid 
in the case of ^^Ar. 



III. PROBING THE CHARGE DENSITY 

Although there are strong indications that the 2si/2 is 
half filled in ^^Ar, we choose this nucleus to illustrate ex- 
perimental signals associated with the bubble occurrence, 
as explained above. It should be noted that similar con- 
clusions could be drawn on ^^Ar or other bubble nuclei 
with respect to the present calculations. ^^Ar is an un- 
stable nucleus which can be produced in present exotic 
beam facilities with a typical intensity of 10'^ pps. Several 
experimental studies have already been performed and 
we will analyze in the following these data with respect 
to a bubble manifestation. The most direct observable 
related to a proton bubble is the charge density. 



A. The charge radius 

Charge radii of Ar isotopes have been precisely mea- 
sured from -^^Ar to ^^Ar i27j,|28|, using laser spectroscopy 
of fast beams at the ISOLDE isotope separator. The ob- 
tained value for ^^Ar is Rc=3.44 ± 0.01 fm. However, the 
calculated "^^Ar r.m.s. charge radii with the HF or HFB 
models are mainly sensitive to the surface part of the 
proton density, and are not adequate quantities to detect 
a bubble in the center of the nucleus: a typical variation 
of 0.02 fm is found between calculations assuming the 
bubble and the no-bubble hypothesis. This difference 
is not enough to lead to a clear signature of a bubble 
occurrence. It is therefore necessary to look for comple- 
mentary observables which could be more sensitive to a 
bubble manifestation. 



B. Electron scattering 

The ideal experiment to probe the charge density 
would be electron scattering. Presently electrons scat- 
tering on unstable nuclei like ^^Ar is not possible, but 
it is expected in next generation facilities such as RIBF 
in Riken \29\ of FAIR in GSI In this case accel- 

erated electrons would scatter on a radioactive beam of 
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FIG. 3: Angular dependence of the form factor obtained for 
300 MeV electron scattering on *® Ar, using either the HF-SkI5 
density (dashed line) or the HF-SLy4 density (solid line) 

^^Ar kept in a storage ring. Such an experiment may be 
feasible in the next decade. It is therefore worthwhile to 
check how accurately the electron scattering experiment 
could probe the proton bubble presence. We have per- 
formed calculation for 300 MeV electron scattering ^^Ar. 
Fig. [3] shows the angular dependence of the form factor 
defined by : 

Fiq) = J Pc{r)e:'-''dr (2) 

where pc is the charge density calculated with HF 
model using SLy4 or SkI5 interactions, and q is the trans- 
ferred momentum, related to the incident momentum p 
and the scattering angle 6: 

e 

9 = 2psm- (3) 

Both the no-bubble and the bubble cases exhibit a 
diffraction pattern. However, above 45 deg. the angular 
distribution is located at higher magnitude in the case of 
a bubble than without a bubble. Around 80 deg. an an- 
gular shift between the two distributions is noticed, the 
minima of the bubble distribution being 10 deg. larger 
than the normal one. At 90 deg. the two angular distri- 
butions are in opposite phases, showing that this angle 
is the optimum one to disentangle between the two hy- 
pothesis. At larger angle the angular shift increases but 
the magnitude becomes too small for a clear measure- 
ment. In summary, an optimal experiment to search for 



a bubble in very neutron-rich Argon isotopes would be 
elastic electron scattering at 300 MeV, focused on scat- 
tered electrons between 40 degrees and 100 degrees. 



IV. DIRECT REACTIONS AND COLLECTIVE 
MODES 

As stated above, such an experiment with electron 
scattering on exotic beams is not feasible yet. It is there- 
fore interesting to analyze what hints direct reactions 
could provide with respect to a bubble manifestation. 
It is well known that direct reactions probe the surface 
of the nucleus, so it is not possible to directly detect a 
bubble located in the center of the nucleus. However, in 
the Argon isotopes, the bubble is a straightforward con- 
sequence of the inversion between the ld3/2 and the 2si/2 
levels, provided that the 2si/2 remains depopulated. Di- 
rect reactions are the most accurate tool to study these 
features. For instance transfer reactions such as (d,^He) 
would allow to measure the spectroscopic factors associ- 
ated to the 2si/2 level. Collective modes such as giant res- 
onances and low-lying states may also provide useful in- 
formation about a bubble manifestation. We shall again 
illustrate the bubble occurrence on ^^Ar, which could 
more probably occur on more neutron- rich isotopes. 



A. Giant resonances 

Giant resonances are collective modes involving all the 
nucleons. To disentangle the effect of the bubble from the 
one due to the interaction, namely due to the modifica- 
tion of the single particle spectrum, HF and HFB calcu- 
lations have also been performed with both the SLy4 and 
the SkI5 interactions in the ''^Ca case, where no bubble 
effect is present. 

In order to perform such a study and compare with 
the data, HFB-I-QRPA calculations have been performed. 
All detailed presentation of the model can be found in 
Ref. [Slj. The residual interaction is derived from the 
Skyrme functional, as well as the pairing residual inter- 
action, derived from the pairing part of the functional 
used in the HFB calculation. We have first checked that 
the HFB+QRPA results using SkI5 give results similar 
to the HF-I-RPA using SLy4, that is the no-bubble hy- 
pothesis. Therefore it is sufficient to perform calculations 
using SLy4 without pairing effects. The Landau-Migdal 
approximation is used, leading to a small breaking of the 
self-consistency, cured by a renormalization of the resid- 
ual interaction to set the center of mass spurious mode 
to zero energy fsij. The typical renormalization values 
are 0.9. 

The main effects on giant resonances of ^^Ar which 
are not present in ^^Ca are obtained in the low energy 
part of the dipole and monopole resonances (Fig [5]). No 
strong discrepancy is observed between the bubble and 
the normal cases in this low energy area. In the case of 
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FIG. 4: Isoscalar monopole and isovector dipole strength 
functions obtained with the HF+RPA approach and the SkI5 
interaction (solid line) or the SLy4 interaction (dashed line) 
in 4«Ar 



the monopole response, the lowest energy state is shifted 
by 1.5 MeV at higher energy in the bubble case, from 
14.5 MeV, to 16 MeV. It may be interpreted as a slight 
increase of the nucleus incompressibility for a soft mode 
in the bubble case: protons are depleted in the center of 
the nucleus and are redistributed in the surface, making 
more difficult a soft compression of the nucleus. 

The isovector El strength exhibits an increase of the 
low- lying part, around 13 MeV, in the bubble case. This 
effect is not present in the SkI5 calculations on ^^Ca: 
the strength enhancement is more important in the "^^Ar 
case than in the "'^Ca case. This soft mode may be in- 
terpreted as the oscillation of neutron with respect to 
protons which are surface localized. 

However, the above mentioned effects are hardly siz- 
able, especially for giant resonances related effects, where 
the spreading width implies a large experimental width 
which can mask these effects. 



B. Low- lying states 

The low- lying states of "^^ Ar have been recently studied 
using direct reactions, either by Coulomb excitation [33 
or by proton scattering in inverse kinematics 



33|. As 



stated above, they are an accurate tool to study the level 
inversion as well as the eventual population of the 2si/2 
state due to the pairing effect. 

The study of the first 2+ state of ^^Ar by proton scat- 
tering shows a significant contribution of the neutrons to 



the excitation [33*1 . Hence this state is not the best probe 
for the proton shell structure. However, the study of the 
proton contribution through the B(E2) value could pro- 
vide relevant information. Fig. O shows the isoscalar 
quadrupole response of ^^Ar, calculated by HF+RPA 
with the SLy4 interaction (no bubble) and the SkI5 one 
(bubble) . There is a strong enhancement of the strength 
of the 2f state in the bubble case. The predicted B(E2) 
values strongly differs: 24 e^.fm'' in the normal case and 
256 e^.fm"* in the bubble case. This last case is in very 
good agreement with the data [H : B(E2)=218 ± 31 
e^.fm"*. This enhancement is partly due to the inversion 
between the ld3/2 and 2si/2 states: more proton particle- 
hole configurations contribute in the bubble case than in 
the normal case. However, the main effect comes from the 
RPA residual interaction, which is very different in the 
two cases for this low-lying mode. Namely, the density- 
dependent terms of the residual interaction could explain 
this change, since the density profile is different in the two 
cases. Moreover as stated above, HFB-I-QRPA calcula- 
tions using the SkI5 interaction, give similar results to 
the normal case, showing that it is the density depen- 
dence of the residual interaction which drives the B(E2) 
enhancement. These terms are generated by the density- 
dependent part of the effective nucleon-nucleon interac- 
tion itself. Therefore the study of bubble nuclei could not 
only provide information on the spin-orbit term as stated 
in Rcf. [l3)[iB| , but also on the importance of the density- 
dependent term of the nucleon-nucleon interaction. The 
measurement of the B(E2) in bubble nuclei may open a 
way to explore the impact of this term in the energy den- 
sity functional. This shows that the manifestation of a 
bubble in ''^Ar might be necessary to reproduce the data, 
namely the first 2+ state. 



V. CONCLUSIONS 

We have investigated the bubble occurrence in nuclei, 
due to an inversion of the Si/2 state with the state usu- 
ally located above. The neutron-rich Argon isotopes are 
the best candidates of the nuclear chart, since the pro- 
ton 2si/2 state largely contributes to the density in the 
center of the nucleus, and because the inversion effect 
is known to increase with the neutron number. Very 
neutron-rich nuclei such as ^*Ar are optimal candidates. 
In less neutron-rich nuclei such as ^^Ar, the situation is 
less favorable since the (2si/2,ld3/2) inversion is reduced, 
implying a significant occupation probability for the 2si/2 
state, leading to a weakening of the bubble effect. 

*^Ar is chosen as an illustrative case, considering both 
the bubble and no-bubble hypothesis: experimental data 
are available on this nucleus, and the predictions are sim- 
ilar for more neutron-rich isotopes. The p-states neu- 
tron splitting may be related to a proton depletion, and 
the B(E2) is also very well reproduced in *^Ar, consid- 
ering the bubble hypothesis. These indications call for 
a more detailed experimental investigation in this nu- 
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cleus, as well as more neutron-rich ones, which could be 
undertaken by elastic electron scattering with the next 
generation exotic beam facilities such as FAIR or RIBF. 

The bubble effect is weakened by pairing effects which 
can populate the 2si/2 state. However the inverted s 
and d level spacing increases for very neutron-rich Argon 
nuclei, and even in the case of a strong pairing, some 
isotopes more neutron-rich than ^^Ar should exhibit a 
bubble, as predicted in the dramatic case of Ar. Hence, 
direct reactions on so-S^^j. next generation radioactive 
beam facilities will also be of great interest. They should 
focus on low-lying states since giant resonances do not 
exhibit strongly different patterns in the case of a bubble, 
compared to the normal case. 

The discovery of bubbles in nuclei is certainly an im- 
portant issue for nuclear structure. This study calls for 
other theoretical investigations, as well as an important 
experimental program devoted to frontline exotic nuclei, 
such as neutron-rich argon isotopes. 



FIG. 5: Isoscalar quadrupole strength function obtained with 
the HF-I-RPA approach and the SkI5 interaction (solid line) 
or the SLy4 interaction (dashed lines) in *^Ar 



Acknowledgments The authors thank F. Leblanc, O. 
Sorlin and D. Verney for fruitful discussions. 



[1] H.A. Wilson, Phys. Rev. 69 (1946) 538. 

[2] P. Siemens, H.A. Bethe, Phys. Rev. Lett. 18 (1967) 704. 

[3] W.J. Swiatecki, Phys. Script. 28 (1983) 349. 

[4] O. Bohigas, X. Campi, H. Krivine, J. Treiner, Phys. Lett. 

B64 (1976) 381. 
[5] J. Decharge, J.-F Berger, M. Girod, K. Dietrich, Nucl. 

Phys. A716 (2003) 55. 
[6] Y. Yu, A. Bulgac, P. Magierski, Phys. Rev. Lett. 84 

(2000) 412. 

[7] M. Bender, K. Rutz, P.-G. Reinhard, J.A. Maruhn, W. 

Greiner, Phys. Rev. C60 (2003) 034304 
[8] J.M. Cavedon et al, Phys. Rev. Lett AQ (1982) 978. 
[9] X. Campi and D.W.L. Sprung, Phys. Lett. B46 (1973) 

291. 

[10] M. Beiner and R.J. Lombard, Phys. Lett. B47 (1973) 
399. 

[11] A. Ozawa, T. Kobayashi, T. Suzuki, K. Yoshida, I. Tani- 

hata, Phys. Rev. Lett 84 (2000) 5493. 
[12] E. Becheva et al, Phys. Rev. Lett 96 (2006) 012501. 
[13] Z. Elekes et al, Phys. Rev. Lett Q8 (2007) 102502. 
[14] B.G. Todd-Rudel, J. Piekarewicz, P.D. Gottle, Phys. Rev. 

C69 (2004) 021301(R). 
[15] J. Piekarewicz, J. Phys. G34 (1996) 467. 
[16] M. Grasso, Z. Ma, E. Khan, J. Margueron, Nguyen Van 

Giai, Phys. Rev. C76 (2007) 044319. 
[17] C.A. Ogilvie et al, Nucl. Phys. A465 (1987) 445 
[18] P.G. Reinhard and H. Flocard, Nucl. Phys. A584 (1995) 

467 

[19] L. Gaudefroy et al, Phys. Rev. Lett. 99 (2007) 099202 
[20] E. Chabanat, P. Bonche, P. Haensel, J. Meyer, R. Scha- 
effer, Nucl. Phys. A635 (1998) 231. 



[21] R.B. Firestone, Ta6Ze o/ isotopes 1996, Wiley Interscience 
Ed. 

[22] T.R. Werner et al., Nucl. Phys. A597 (1996) 327 

[23] G.A. Lalazissis, D. Vretenar, P. Ring, M. Stoitsov, and 

L. M. Robledo, Phys. Rev. C60 (1999) 014310 
[24] S. Hilaire, M. Girod, http://www-phynu.cea. fr/^ 

science_enJigne/carte_potentiels_microscopiques/ 

carte_potentieLnucleaire.htm (2006) 
[25] R. Rodriguez-Guzman, J.L. Egido, L.M. Robledo, Phys. 

Rev. C65 (2002) 024304 
[26] L. Gaudefroy et al, Phys. Rev. Lett. 97 (2006) 092501. 
[27] A. Klein et al., Nucl. Phys. A607 (1996) 1 
[28] K. Blaum et al., Hyperfine Interactions 162 (2005) 101 
[29] T. Motobayashi, Recent Achievements and Perscpetive 

in Nuclear Physics, Proc. of the Vth Italy- Japan Sympo- 
sium, Naples (2005) 563 
[30] B. Rubio and T. Nilsson, Nucl. Phys. News 16 (2006) 5 
[31] E. Khan, N. Sandulescu, M. Grasso, N. V. Giai, Phys. 

Rev C66 (2002) 024309. 
[32] A. Gade et al., Phys. Rev. C68 (2003) 014302 
[33] L.A. Rileyet al, Phys. Rev. C72 (2005) 024311 



